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Abstract. - We numerically investigate the nature of the phase transition of the XY model in the 
heptagonal lattice with the negative curvature, in comparison to other interaction structures such 
as a flat two-dimensional (2D) square lattice and a small-world network. Although the heptagonal 
lattice has a very short characteristic path length like the small-world network structure, it is 
revealed via calculation of the Binder's cumulant that the former exhibits a zero-temperature 
phase transition while the latter has the finite-temperature transition of the mean-field nature. 
Through the computation of the vortex density as well as the correlation function in the low- 
temperature approximation, we show that the absence of the phase transition originates from the 
strong spinwave-type fluctuation, which is discussed in relation to the usual 2D XY model. 



A common belief in the complex network researches is 
that critical phenomena in networks often exhibit mean- 
field natures. For the small-world network structure pro- 
posed by Watts and Strogatz (WS) [1], it is well-known 
that statistical mechanical model systems such as the Ising 
and the XY models have phase transitions completely dif- 
ferent from their counterparts in one- and two-dimensional 
(2D) regular lattices [2,3]. Furthermore, the universality 
class of the phase transition has been known to undergo 
a dramatic change as soon as the density of shortcuts has 
a finite nonzero value, which is accompanied by the struc- 
tural phase transition from the large world to the small 
world. These existing findings imply the importance of 
the topological interaction structure in the determination 
of the universality class of the phase transition. Recently, 
the Ising model in the heptagonal lattice, which has com- 
plicated geometry with a negative curvature [4] , has been 
shown to possess the mean-field critical exponents, consis- 
tently with the common belief [5,6]. 

In the present Letter, we take the one of the mostly 
well-studied statistical mechanical model systems, the XY 
model, and study the nature of the phase transition in the 
heptagonal lattice. The comparison with the usual flat 
2D square lattice sheds a light on the role of the geome- 
try with different curvatures, while the comparison with 
the small- world network structure indicates that, different 
from a naive expectation, the short path length cannot 




Fig. 1: Schematic representation of a heptagonal lattice upto 
the 4th level (L = 4), projected on the Poincare disk [7]. The 
surrounding circle corresponds to the points at infinity, and 
the lattice within it is composed of congruent heptagons with 
respect to the negative curvature metric [8] . 



solely determine the universality class. Extensive Monte- 
Carlo computations of the Binder's cumulant and the vor- 
tex density, combined with the numerical calculation of 
the spin-spin correlation function, bring us to the con- 
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Fig. 3: (Color online) Ferromagnetic order parameter m of 
heptagonal lattices, as a function of the temperature T in units 
of J/ks- The spins in the two outmost layers are excluded in 
computing m to remove the undesirable boundary effects. 



Fig. 2: (Color online) The characteristic path length t versus 
the system size N for the 2D square lattice, heptagonal lattice, 
and the WS network (from top to bottom). Both the heptag- 
onal lattice and the WS network exhibit I ~ In AT, while the 
2D square lattice shows I ~ N 1 ^ 2 . The latter is more clearly 
seen in the inset in which the same data points are plotted in 
log-log scales. 



elusion that the absence of the finite-temperature phase 
transition in the heptagonal lattice originates not from 
the vortex fluctuation, but from the strong spinwave fluc- 
tuation. 

We start from a brief description of the heptagonal lat- 
tice [5,6] (see Fig. 1). Suppose that we start construct- 
ing the heptagonal lattice from a single heptagon (the 1st 
level), the central one in Fig. 1. Its seven nearest neighbor 
heptagons constitute the 2nd level, and the next nearest 
neighbors the 3rd level, and so on. We denote the set of 
heptagons in the Zth level as H(l) and call set of vertices on 
the outward boundary of H(l) as the lih layer. The total 
number N(L) of spins in the heptagonal lattice containing 
all the layers upto Lth one has been shown to increase ex- 
ponentially with L [5,6], which implies the following two 
interesting features: (i) The surface- volume ratio does not 
vanish but remains finite in the thermodynamic limit, (ii) 
The average path length £ increases logarithmically with 
the system size N, i.e., £ ~ In N, since the spins separated 
by a long distance on a given layer are connected by much 
shorter path through inter-layer connections. In Fig. 2 we 
plot I versus N for the heptagonal lattice, the 2D square 
lattice, and the WS small- world network. It is clearly seen 
that £ ~ IniV for the heptagonal lattice and for the WS 
network, and £ ~ y/N for the 2D square lattice. From 
the above comparisons, one may naively expect the same 
universality class both for the heptagonal lattice and the 
WS network, however, it is clearly shown below that this 
is not the case. 

The XY model is one of the most important model sys- 
tems not only in statistical physics but also in condensed- 
matter physics due to the existences of real-world exam- 



ples, and is described by the Hamiltonian [9]: 

H = -J^2 cos(0j - <f>j), (1) 
(v) 

where J is the coupling constant, the sum is over all near- 
est pairs in the system, and the variable <f>i at the ith site 
can represent either the direction of the two-dimensional 
spin in magnetic systems or the phase of the complex 
order parameter in superconducting systems [10]. In a 
two-dimensional lattice, the Mermin- Wagner theorem [11] 
states that there cannot be a true long-range order at fi- 
nite temperatures, however, it has been known that the 2D 
XY model exhibit a quasi-long-range order at low temper- 
atures, characterized by the algebraic decay of the corre- 
lation function. 

We perform Monte Carlo simulations based on the 
Metropolis algorithm. The ferromagnetic order parameter 
m = (1/N) J2j e 1 is computed and the Binder's cumu- 
lant for the system of the size N 



is measured with the ensemble average (• • •). We plot (|m|) 
in Fig. 3 as a function of the temperature T in units of 
J /kg with the Boltzmann constant kg, which clearly in- 
dicates that the heptagonal lattice lacks a true long-range 
order, qualitatively similarly to the 2D regular square lat- 
tice. The lowest excitation in the XY model is the gapless 
spinwave excitations, and thus there is no essential singu- 
larity at zero temperature in thermodynamic quantities 
like susceptibility. The next question for the presence of a 
quasi-long-range order can be answered by observing the 
Binder's cumulant, Un- Figure 4 shows Un versus the 
temperature T for (a) the heptagonal lattice, (b) the 2D 
square lattice, and (c) the WS small-world network. As 
were already found in the literature, the 2D square lattice 
shows merging of Un for different TV's in the whole low- 
temperature phase with the quasi-long-range order [12], 
while the WS network exhibits a unique crossing at a 
well-defined critical temperature splitting the truly or- 
dered phase and disordered phase [3] . The merging of the 
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Fig. 4: (Color online) Binder's cumulant Un versus the tem- 
perature T for (a) heptagonal lattices, (b) 2D regular square 
lattices, and (c) WS networks (Data in Ref. [3] are used). In 
(a), we have excluded two outmost layers in the computation 
of Un to avoid surface effects. Inset of (a): expansion of the 
temperature range T < 0.1, showing that Un's neither cross 
nor merge. 



Binder's cumulant for the 2D square lattice is due to the 
diverging correlation length in the whole low-temperature 
phase, which has also been detected by the merging of the 
ratios of two correlation functions [13]. 

For the heptagonal lattices [see Fig. 4(a)], we exclude 
spins in outmost two layers in the computation of m for 
a given system size in order to avoid the artifact origi- 
nating from the nonvanishing surface-volume ratio in the 
thermodynamic limit. The Binder's cumulant for heptag- 
onal lattices in Fig. 4(a) displays a completely different 
behavior than for 2D regular lattice and WS networks: 
The curves for Un obtained from different sizes do not 
merge even at very low temperatures T < 0.1 [see Inset 
of Fig. 4(a)]. Consequently, we conclude that the XY 
model in the heptagonal lattice is disordered at any finite 
non-zero temperature and does not even possess the quasi- 
long-range order present in the 2D regular lattice. This is 
a remarkably surprising result from the viewpoint of the 
topological interaction structure: In contrast to the WS 




Fig. 5: (Color online) Vortex number density n v in a heptag- 
onal lattice with seven layers as a function of T is measured 
for each level I — 1, 2, • ■ ■ , 7. The proliferation of vortices oc- 
curs around T = 0.4, indicating that the absence of order in 
Fig. 4(a) is not from the vortex fluctuation. 



network, the short path length in the heptagonal lattice 
does not increase the effective dimension of the system 
making the ordered state more probable. We next pursue 
the answer to the question on the origin of the enhanced 
fluctuation in the heptagonal lattice. 

The quasi-long-range order present in the 2D XY model 
is known to be broken by the proliferation of free vor- 
tices [9], Accordingly, the above found absence of the 
quasi-long-range order in the heptagonal XY model at 
finite temperatures could be explained by the strong vor- 
tex fluctuation even at extremely low temperatures. In 
order to test this possible scenario, we measure the vortex 
number density n v (l) for the Ith level heptagons defined 
as 



n v (l) = 




(3) 



peH(i) 



where |-ff(0l is the number of heptagons in the Ith level, 
the summation ^ p is taken in the counter-clockwise direc- 
tion over the bonds surrounding the heptagon p, and the 
phase difference 4>i — <j)j is defined modulo 2tt. Although n v 
cannot detect the difference between the free and bound 
vortices, it reflects the Kosterlitz-Thouless (KT) transition 
since it changes from zero to nonzero smoothly around the 
transition temperature in the 2D XY model. We in Fig. 5 
plot n v as a function of T for the heptagonal XY model 
composed of seven layers (L = 7): Vortices are observed to 
be generated at some finite temperature around T = 0.4 
and the number density for each layer does not show a 
significant difference. Overall, one sees that the tempera- 
ture region around T = 0.4 in which thermal vortices are 
generated is well above the zero temperature, which leads 
us to the conclusion that the absence of any long-range 
order (be it genuine or quasi) found in the heptagonal 
XY model cannot be explained from the vortex degree of 
freedom. Accordingly, we below focus on the spinwave ex- 
citation in the XY model in the low-temperature regime. 

The spin-spin correlation function within the low- 
temperature approximation, i.e., cos 6* « 1 — 2 /2, is easily 
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Fig. 6: (Color online) The low-temperature correlation func- 
tion is shown in the form of g(d), which is defined as the aver- 
age of gjk for all pairs of (j, fc)'s separated by the given shortest 
path length d [see Eqs. (4) and (5)]. In (a), g increases linearly 
with d for the heptagonal lattice (with L = 6), indicating the 
exponential decay of the correlation function. In (b) the same 
data are plotted in lin-log scales, and g increases logarithmi- 
cally for the 2D square lattice (N — 40 x 40). It is seen that 
for the WS network (TV = 1600) g(d) saturates toward a finite 
value as d is increased, indicating the existence of the ordered 
phase at low temperatures. 



computed [14] to yield 



C jk 
9jk 



(1/2) (G 



:)\ 



s e -9ik/T 

Gjk + Gkk 



G 



kj i 



(4) 
(5) 



where the lattice Greens function Gjk is the inverse of the 
discrete Laplacian [15] defined by Gjf} = kjSjk — A/fc with 
the degree kj and the adjacency matrix Ajk [16]. It is 
natural to measure the distance between the two spins at 
the jth and the fcth vertices by the shortest path length, 
which is especially important for interaction structures not 
put on geography, like the WS network. If there exist 
paths of different lengths, we believe that the path with 
the shortest length must give the most contribution to the 
correlation function. 

In Fig. 6, we present spatial average g(d) of gjk in Eq. (5) 
taken for all pairs split by the shortest path length d. The 
logarithmic increase of g(d), and accordingly the algebraic 
decay of C(d) in Eq. (4), for the 2D square lattice clearly 
shown in Fig. 6(b) indicates the existence of the quasi- 
long-range order in the low-temperature phase below the 
KT transition [9] . In Fig. 6 it should be noted that the ab- 
sence of the long-range order is reflected as the divergence 
of g(d) as d is increased. For the WS network, on the 
other hand, g{d) does not diverge, but saturates to a well- 



defined value, which is in an agreement with the expec- 
tation from Ref . [3] : The existence of the true long-range 
order of the XY model in the WS network at sufficiently 
low temperatures should be reflected as a finite and non- 
vanishing correlation function Cjk as djk is increased. In 
contrast, the behavior of g(d) for the heptagonal lattice in 
Fig. 6(a) is completely different both from the 2D square 
lattice and the WS network: g(d) increases with d like the 
2D square lattice [shown in Fig. 6(b)], however, not log- 
arithmically but linearly. The linear increase of g with d 
implies that the correlation function C (d) decays exponen- 
tially, which clearly implies the absence of any long-range 
order at sufficiently low temperatures. As a consequence, 
we conclude that the absence of the ordered phase in the 
heptagonal XY model originates not from the vortex fluc- 
tuation but from the spinwave excitation, which is strong 
enough to destroy even the quasi-long-range order present 
in the low-temperature KT phase of the 2D square XY 
model. 

We emphasize that our observation of the absence of 
the ordered phase in the heptagonal XY model is con- 
trasted to the report for the heptagonal Ising model [5,6]: 
The geometrical difference between the planar 2D lattice 
and the heptagonal lattice with a negative curvature only 
changes the critical exponents of the Ising model to the 
mean-field ones. In other words, for the Ising model the 
overall strength of fluctuation appears to be reduced by 
an increase of an effective dimensionality, reflected as the 
small- world behavior of very short path lengths, similarly 
to the WS networks. In a sharp contrast, the critical phe- 
nomena of the XY model at low temperatures depend 
sensitively on how the correlation function behaves, which 
appears to lead completely different critical behaviors. It 
is also interesting to note that the very same fact that all 
the vertices are closely connected by short path lengths 
yield very different behaviors, i.e., stimulating the mean- 
field nature for the Ising model, and the enhancement of 
spinwave fluctuation to remove any ordered phase in the 
XY model. This is not only an example that the effect 
of lattice geometry can solely change the criticality to a 
great extent, but also that the close connectedness does 
not guarantee the mean-field criticality. 

We believe that the absence of an ordered phase in the 
negatively curved system should apply in a broad range 
of different cases, since it is based on the quadratic ap- 
proximation of the interaction which must be valid for a 
variety of different interaction forms, including Villain for- 
mulation of the XY model. Furthermore, even when the 
quadratic approximation fails, we believe the absence of 
an order should still be true, since it appears to depend 
strongly on the topological lattice structure, rather than 
the detailed form of the interaction: The exponentially 
decaying correlation function must be related to the fact 
that the boundary increases exponentially with the dis- 
tance from any point. Since it is the basic geometric fea- 
ture of a negatively curved surface, not restricted to the 
heptagonal structure, the physics will remain unaltered 
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also in the continuum. 

In summary, we have investigated the XY model in hep- 
tagonal lattices in parallel to studies of 2D square lattices 
and the WS small-world networks. Extensive Monte-Carlo 
simulations have been performed to compute the Binder's 
cumulant, which has shown clearly that neither the true 
long-range order nor the quasi-long-range KT order is es- 
tablished at any finite temperatures. The vortex number 
density has been shown to be not the origin of the ab- 
sence of an ordered phase. From the numerical calcula- 
tion of the correlation function we have found that the 
short path lengths facilitate the strong spinwave fluctua- 
tion, which eventually removes any order, leading to the 
exponentially decaying correlation function at any nonzero 
temperature. We have also discussed the recent works on 
the Ising model in the same heptagonal lattice. 
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